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Through an extensive numerical study, we find that the low-lying, quasi-degenerate eigenenergies 
of weakly-interacting trapped A'^ bosons with total angular momentum L are given in case of small 
L/N and sufficiently small L hy E = Lhu -\- g[N{N - L/2 - 1) -\- 1.59 n(n - l)/2], where uj is the 
frequency of the trapping potential and g is the strength of the repulsive contact interaction. The 
last term arises from the pairwise repulsive interaction among n octupole excitations and describes 
the lowest-lying excitation spectra from the Yrast line. In this case, the quadrupole modes do not 
interact with themselves and, together with the octupole modes, exhaust the low-lying spectra that 
are separated from others by A-linear energy gaps. 



The Yrast state, which is a subject of active research in 
nuclear physics , is the lowest-energy state of a system 
of particles for a given total angular momentum (AM). 
This state is important in that the whole of the excita- 
tion energy is used up for the rotation of the system and 
hence the system is at zero temperature, opening up the 
possibility of performing precise spectroscopic measure- 
ments of energy levels close to the Yrast state. Recently, 
Mottelson |^ has pointed out that similar problems arise 
in rotating Bose-Einstein condensates (BECs) of trapped 
atomic vapor Isl . Two recent observations of vortices in 
trapped BECs |H| have lent agreat impetus to theoretical 
studies of this subject 

The Yrast state is obtained by distributing a given 
AM L over N bosons so as to minimize the total en- 
ergy. When the trapping potential is harmonic, there 
would be a huge number of such partitions for L ^ 1 
and hence a huge degeneracy, were it not for the inter- 
actions The problem thus reduces to finding how a 
weak interaction lifts the degeneracy, selects the Yrast 
state, and determines the low-lying excitations from this 
state. Wilkin et al. Q have found that when the interac- 
tion is attractive, all the AM in the Yrast state resides in 
the center-of-mass motion of the system. Mottelson ||] 
discussed that low-lying excitations are well described by 
collective modes whose excitation operators are given by 




(1) 



where A is an integer, and Xp and j/p are the position op- 
erators of the p-th particle. The corresponding excitation 
energy is given by 
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where g denotes the strength of the contact interaction. 
Equation (^ suggests that for the case of attractive inter- 
actions (i.e., g < 0), the AM of the Yrast state is carried 



by the dipolar (A = 1) mode, in agreement with Ref. 
while for the case of repulsive interactions it is carried by 
the quadrupole (A — 2) or octupole (A = 3) modes [D, 
as they have the greatest energy gain per unit of AM. A 
more elaborate study taking the mode-mode interaction 
into account has shown that the quadrupole modes have 
slightly larger energy gains than the octupole modes [Q. 

Bertsch and Papenbrock |^ have performed numeri- 
cal diagonalization for small systems {N — 25 and 50), 
finding that the energy of the Yrast line for L < iV can 
be given by i? = Lhuj + gN{N — L/2 — 1), where cu is 
the frequency of the trapping potential. The correspond- 
ing eigenstate was very recently shown to exist (see 
also Refs. [^^). This state partitions the AM equally 
among the particles j|,|ri|,|l^ and smoothly crosses over 
to the many-body single vortex state |^ as i approaches 
N. 

In this paper, we report on the results of our extensive 
numerical study of this system with up to 640000 parti- 
cles and for a total AM up to 30. Our primary finding is 
that of the lowest-lying excitation spectra from the Yrast 
line that arise from the pairwise repulsive interaction be- 
tween octupole modes; the corresponding interaction en- 
ergy is given by 1.59 gn(n — l)/2, where n is the number 
of excited octupole modes . These energy levels are 
separated from the other excitations, including the one 
discussed in Ref. Q, by iV-linear energy gaps. We have 
also found that for small L/N and sufficiently small L all 
the many-body eigenstates have surprisingly large over- 
lap (~ 0.99) with trial wavefunctions defined by 
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with Xnx ~ L, 
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where n\ is a non-negative integer and |0) describes the 
exact many-body ground state of the system with L = 0. 

The model we study is the same as that of Refs. [|-||. 
We consider a system of weakly-interacting bosons sub- 
ject to a given AM and trapped in a parabolic confining 
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potential. The system is isotropic in the x — y plane 
and strongly confined in the z direction so that all the 
particles occupy the lowest-energy state in this direction. 
The problem is thus essentially two-dimensional. The 
many-body Hamiltonian for this system can be written 
a.s H — J2i + Here hi denotes the single-particle 
Hamiltonian for the i-th particle 
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where M is the mass of the boson, lo is the frequency 
of the trap, and V = (inh^g/M) '£t<j " ^j) de- 

scribes the interaction between particles {g is the dimen- 
sionless 'scattering length'). In this paper we consider 
the case of g > 0, i.e., repulsive interactions. 

The single-particle spectrum is given by En^m = 
hu {2n + \m\ + 1), where n and m denote the radial 
and AM quantum numbers, respectively. When only low- 
lying states of the many-boson system are concerned, we 
may assume that all the particles are in states with n — 0. 
The single-particle states are then labeled only by m and 
described by (j)m{z) = (z™/\/ 7rm!) exp(— |zp/2), where 
z = X + iy, and the lengths are measured in units of 
(n/Afw)i/2. We use these single-particle states as a basis 
set to rewrite the many-body Hamiltonian as 
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where g = ghuj, the operators and bl^ annihilate and 
create one boson in the single-particle state (j)m, respec- 
tively, and V,nim2m3mi IS the two-body matrix element 
given as |9|] 
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Given a total AM i, the state may be spanned by 
the Fock states |rio, ?^L) with X]i=o l ^« ~ ^ 



and Ej=o,lJ'^ 



L, where iii denotes the occupa- 



tion number of the i-th single particle state 4>i, i.e., 
b\bi\no,ni, . . . ^ul) = nj|no, ni, . . . , ul). Numerical di- 
agonalization of the two-body interaction V is performed 
within the Hilbert subspace subject to these constraints. 

The energy eigenvalue corresponding to the i-th many- 
body eigenstate with total AM L may be written as 
LfiLo -t- eL.i, where cla denotes the interaction energy. 
Because the lowest-energy eigenvalue with total AM L 
{2 < L < N) is given by el = gN{N - L/2 - 1), 
the unique state with L — which is composed of 
A'' bosons in a single-particle state is the ground 
state for gN <C Tilo, and it will therefore be denoted 
as |G). Our trial wavcfunctions are constructed by act- 
ing on the ground state |0) with the collective opera- 
tors, {Q\\, whose second quantized forms are given by 
Qa = (1/V]V)E„. V(™ + A)!/m!A!6j,,+,fe„. 
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FIG. 1. The excitation spectra for iV = 10000, L < 8. The 
energy is measured in units of gN, and the Yrast line is shown 
as a solid line. The collective excitations, Qa|0) (A < 8), are 
linked by a dashed curve, and their energy values are given 
in Table I. Except for the state Qi|0), the states that in- 
clude excitations of the center-of-mass motion are omitted. 
Each state has a very large overlap (> 0.994) with the cor- 
responding trial wavefunction. These states are given in the 
ascending order of energy as Qi|0) for L = 1; Q2IO) for L — 2; 
QslO) for L = 3; (Q2)^|0>, Q^O) for L = 4; QaQalO), QslO) 
(Q2)'|0), (Q3)'|0), Q2Q4IO), QslO) for L = 6; 



for L = 5; 

(Q2)'Q3|0), 
O2(Q3)'|0>, 

for L = 8. 

[(Q2)''|0), (32(Q3)^|0)] for L = 8, are quasi-degenerate, and 
these states are shown on or extremely near the Yrast line. 



Q3Q4IO), Q2Qb\0), QtIO) for L = 7; (Q2)*|0), 
Q3Q5IO), Q2g6|0>, QslO) 
(Q3)^|0)] for L = 6 and 



(Q2)'Q4|0), (Q4)10>, 

The sets, [(Q2)^|0), 



Figure ^ shows the energy spectra for N = 10000 and 
L < 8, where the energy is measured in units of gN. Our 
numerical results show that all the energy eigenstates for 
L < 8 have remarkably large overlaps (> 0.994) with 
the corresponding trial wavcfunctions given in Eq. (|^). 
Among the eigenstates, only "intrinsic" ones, i.e., those 
that do not include dipole excitations, are shown in Fig. |l| 
(the only exception being the state Qi|0) for L = 1), 
as the dipole mode corresponds to the center-of-mass 
motion that does not affect the interaction in a har- 
monic trap. For example, for L = 6 only the "intrin- 
sic" eigenstates, (Q2)^|0), (Q3)^|0), Q2Q4\0), and QslO), 
are shown among all the eigenstates for L = 6, and 
they are identified in the ascending order of energy as 
(Q2)^|0), (Q3)'|0), Q2Q4IO), Q1Q2Q3IO), (Qi)2(Q2)'|0), 
QelO), Q1Q5IO), (Qi)2Q4|0), (Qi)3Q3|0), (Qi)^Q2|0), 
(Qi)^|0). We should note here that two states, {Q2)^\0) 
and ((53)^10)1 are seen as 'a single circle' in Fig. |l] due 
to their quasi- degeneracy. The energy difference between 
them is caused by mode-mode interactions and is on the 
order of g (|0,|l|]. Two states (O2)''|0) and O2(Q3)^|0) 
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for L = 8 are also quasi-degenerate. The Yrast line is 
shown by a solid line in Fig. and the other states are 
separated from this line by energy gaps roughly equal to 
or greater than gN/A. The energy difference gN/A cor- 
responds to the replacement of the factor (Q2)^ included 
in the Yrast states with Q4. 

In Fig.0, the collective excitations discussed by Mottel- 
son, Q\\0) (A = 1, 2, ... , 8), are linked by a dashed curve. 
The energies of these states are given in Table |. These 
values are very well described by Eq. (^; in particular, 
for A < 3 the agreement is perfect within the machine 
accuracy Q. 

For higher (but not very large) AM, the low-lying 
states are again well described by the collective excita- 
tions. In particular, the quadrupole (A = 2) and oc- 
tupole (A — 3) modes well describe the low-lying, quasi- 
degenerate states as shown in Table These quasi- 
degenerate states for 2 < L < 30 are shown in Fig. ||, 
where the energies are measured from the Yrast line and 
are on the order of g. We note that the other excitations 
are separated from these in Fig. 2 by iV-linear energy 
gaps {>gN/A). 

However, the overlap integrals between the collective 
trial wavefunctions and the numerically-obtained low- 
lying eigenstates become small even for small L/N as 
L increases. In fact, in Table [II. the deviations of these 
overlaps from unity seem to include not only a factor of 
L/N but also another factor like L" (s ~ 2) for small 
L. Thus these trial wavefunctions describe the low- lying 
states well for small L/N and sufficiently small L. 

Moreover, for small L/N and sufficiently small L, the 
energies of these quasi-degenerate states are found to be 
described by 1.59 5 x n(n — l)/2, where n is the number 
of excited octupole modes in each state [^Jl4|. This en- 
ergy 'quantization' can also be seen in the data for N = 
640000, L = 18 in Table |l^ as 1.590 ~ 1.59 x 2(2- l)/2, 
9.530 ~ 1.59 X 4(4 - l)/2 and 23.82 ~ 1.59 x 6(6 - l)/2. 
Here, 2, 4 and 6 are the number of excited octupole modes 
in the first, second, and third excited states for L = 18, 
as seen m Table ||. 

These results show that the quadrupole modes do not 
interact with themselves or with the octupole mode in the 
limit oi L/N — > 0, while the octupole modes undergo a 
pairwise repulsive interaction. Althoiigh the mode-mode 
interactions were discussed in Refs. [^,0,|l2| , such concrete 
quasiparticle features of these modes were not pointed 
out previously. Because of such quasiparticle features, 
the Yrast states are well described by ((52)" |0) for even 
L and by (Q2)" QajO) for odd L in that limit (n, n': in- 
tegers). The lowest energy for sufficiently small AM L 
is therefore simply given by the sum of the excitation 
energies of the quadrupole modes (and the energy of an 
octupole mode for odd L). 
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FIG. 2. The low-lying, quasi-degenerate eigenstates for 
A'' = 10000, 2 < L < 30 are shown by open circles. These 
states are well described by the quadrupole and octupole 
excitation modes, as seen in Table II. These energy values, 
which are measured from the Yrast line, are well described 
by 1.59 X n(n — l)/2 in units of g. Here, n is the number 
of excited octupole modes. As the total AM increases, the 
interconversion between these two modes becomes apparent 
(see Table V), which is also suggested by the level repulsion 
between the lowest and first-excited states for larger L. 

As the total AM increases, however, the interconver- 
sion between these two modes becomes significant. This 
is seen in Table ^ and is also suggested by the level repul- 
sion between the lowest and first-excited states for larger 
L in Fig. Although contributions by the modes, Q4, 
Qs and Qe to the Yrast states have been claimed 0, 
the interconversion between modes Q2 and Qs makes a 
greater contribution than these three modes, as seen in 
Table 

Whereas the Yrast state is analytically given for L < N 
[|ll|,|2|, its relations to collective excitations for small 
L /N and sufficiently small L (especially the interconver- 
sion effects between the quadrupole and octupole modes) 
remain unclear. An understanding of these relations may 
help us clarify the integrability in a two-dimensional bose 
system with a harmonic potential that includes weak re- 
pulsive delta- function interactions or a hidden symmetry 
of the Hamiltonian [ p^ . 
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ering the energy contribution by four diagrams in Fig. 5 
in Ref. [Q] up to the order of g. 
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fA/ffiV 





-1 


-3/2 


-1.7499 


-1.8749 


-1.9373 


-1.9686 


-1.9842 



TABLE I. Energies of the collective modes that are excited by acting on the ground state of A'' = 10000 bosons with the 
excitation operator Qx. The agreement with eq. (2) is perfect (within the machine accuracy) for A < 3 and excellent for A > 4. 





lowest 


Ist-excited 


2nd-excited 


3rd-excited 


4th-excited 


5th-cxcited 


L = 9 


0.9951 

(O2)'Q3!0> 


0.9977 

(Q3)^|o> 










12 


0.9804 


0.9815 

(Q2)^g3)'|o) 


0.9959 

(Q3)*|o) 








15 


0.9820 

(Q2)'Q3|0) 


0.9846 

(02)^(Q3)^|0) 


0.9936 

(Q3)^|0) 








18 


0.9374 

(Q2)^|o> 


0.9363 

(Q2)'^(Q3)'|0) 


0.9812 


0.9906 






30 


0.7674 


0.7447 

(Q2)^'(Q3)^|0) 


0.9105 

(Q2)^(Q3)*jO> 


0.9418 


0.9593 

(O2)^(Q3)**|0) 


0.9735 

(Q3)^°|o> 



TABLE II. Values of the overlap integrals between the numerically found low-lying quasi-degenerate states with = 10000 
for various values of L and the corresponding normalized trial wavefunctions (shown on the second line in each cell). These 
values are close to unity for small L/N and sufficiently small L. 
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L = 6 


8 


10 


12 


14 


16 


18 


N = 5000 


0.0049 


0.0117 


0.0227 


0.0385 


0.0596 


0.0861 


0.1179 


10000 


0.0025 


0.0059 


0.0115 


0.0196 


0.0307 


0.0450 


0.0626 


20000 


0.0012 


0.0030 


0.0058 


0.0099 


0.0156 


0.0230 


0.0323 



TABLE III. The deviations of overlap integrals from unity are shown for N = 5000, 10000, 20000 and 
L = 6,8, 10, 12, 14, 16, 18, where each overlap integral is the one between a normalized trial wavefunction with n quadrupole 
excitations {i.e., oc {Q2)"\0) ) and a numerically-obtained lowest-energy state for L = 2n (n: integer). These values seem to 
include not only a factor of L/N but also another factor like L" (s ~ 2) for small L. 





N = 10000 


40000 


160000 


640000 


Ist-excited 


1.721 


1.622 


1.597 


1.590 


2nd-excited 


9.564 


9.538 


9.532 


9.530 


3rd-excited 


23.82 


23.82 


23.82 


23.82 



TABLE IV. Energies (in units of g) of the low-lying quasi-degenerate states with L = 18 for various values of TV. They are 
measured from the Yrast line. As the ratio L/N becomes small, the convergence of these energies can be seen. 





(Q2)^^|o> 


(Q2)^^(Q3)^|0> 




((32)'='Q4|0> 


((32)"'(Q3)'Q4|0> 


lowest 


0.7674 


0.5307 


0.1045 


0.1338 


0.0737 


Ist-excited 


0.5400 


0.7447 


0.1928 


0.0937 


0.1035 



TABLE V. Values of the overlap integrals between the lowest/first-excited states and some trial wavefunctions for N = 10000 
and L = 30. The interconversion between the quadrupole and octupole modes makes a much larger contribution to the Yrast 
state than other modes such as Q4. 
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